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'oq" ■ The masses of the ground state and excited heavy tetraquarks with hidden charm 

and bottom are calculated within the relativistic diquark-antidiquark picture. The 
\ dynamics of the light quark in a heavy-light diquark is treated completely relativis- 

tically. The diquark structure is taken into account by calculating the diquark-gluon 
form factor. New experimental data on charmonium-like states above the open 
charm threshold are discussed. The obtained results indicate that X(3872), Y(4260), 
Y(4360), Z(4433) and Y~(4660) can be tetraquark states with hidden charm. 
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Recently the significant experimental progress has been achieved in heavy hadron spec- 
troscopy. Several new charmonium-like states, such as X(3872), Y (4260), F(4360), Y(4660), 
Qh| Z(4430), etc., were observed [l[ which cannot be simply accommodated in the quark- 
antiquark (qq) picture. These states can be considered as indications of the possible ex- 
istence of exotic multiquark states 0, 0]. Here we briefly review our recent results for the 
masses of heavy tetraquarks in the framework of the relativistic quark model based on the 
quasipotential approach in quantum chromodynamics. We use the diquark-antidiquark ap- 
proximation to reduce a complicated relativistic four-body problem to the subsequent more 
simple two-body problems. The first step consists in the calculation of the masses, wave 
functions and form factors of the diquarks, composed from light and heavy quarks. At the 
final step, a heavy tetraquark is considered to be a bound diquark-antidiquark system. It is 
important to emphasize that we do not consider a diquark as a point particle but explicitly 
take into account its structure by calculating the form factor of the diquark-gluon interaction 
in terms of the diquark wave functions. 

In the quasipotential approach the two-particle bound state with the mass M and masses 
of the constituents mi^ in momentum representation is described by the wave function ^(p) 
satisfying the quasipotential equation of the Schrodinger type 

fb 2 {M) p 2 \ r d 3 q 



*4T(p) = / 7^^V(P> * M )^(q), (1) 
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where the relativistic reduced mass is 

M 4 - (m 2 - ml) 2 

and the on-mass-shell relative momentum squared 

[M 2 - (mi + m 2 ) 2 ][M 2 - mi — 7712 J 



b 2 (M) 



4M 2 



The subscript d refers to the diquark and T refers to the tetraquark composed of a diquark 
and antidiquark. The explicit expressions for the corresponding quasipotentials V^t(p, q; M) 
can be found in Ref. 0]. 
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TABLE I: Masses M and form factor parameters of heavy-light diquarks. S and A denote scalar 
and axial vector diquarks antisymmetric [Q, q] and symmetric {Q, q} in flavour, respectively. 



Quark 
content 


Diquark 
type 




Q = c 






Q = b 




M (MeV) 


£ (GeV) 


C (GeV 2 ) 


M (MeV) 


t (GeV) 


C (GeV 2 ) 


[Q,q] 


S 


1973 


2.55 


0.63 


5359 


6.10 


0.55 


{Q,q} 


A 


2036 


2.51 


0.45 


5381 


6.05 


0.35 


[Q,s] 


S 


2091 


2.15 


1.05 


5462 


5.70 


0.35 


{Q,s} 


A 


2158 


2.12 


0.99 


5482 


5.65 


0.27 



At the first step, we calculate the masses and form factors of the light and heavy diquarks. 
As it is well known, the light quarks are highly relativistic, which makes the v/c expansion 
inapplicable and thus, a completely relativistic treatment of the light quark dynamics is 
required. To achieve this goal we closely follow our consideration of the mass spectra of 
light mesons and adopt the same procedure to make the relativistic potential local by re- 



placing e 1)2 (p) = yJ7ii,2 + P 2 - * •£'1,2 — (M 2 — m 2 1 + m\ 2 )/2M. Solving numerically the 
quasipotential equation (0Q) with the complete relativistic potential, which depends on the 
diquark mass in a complicated highly nonlinear way we get the diquark masses and 
wave functions. In order to determine the diquark interaction with the gluon field, which 
takes into account the diquark structure, we calculate the corresponding matrix element of 
the quark current between diquark states. Such calculation leads to the emergence of the 
form factor F(r) entering the vertex of the diquark-gluon interaction (EJ. This form factor 
is expressed through the overlap integral of the diquark wave functions. Our estimates show 
that it can be approximated with a high accuracy by the expression 

F(r) = 1 - e -f-fr 2 . (2) 

The values of the masses and parameters £ and £ for heavy-light scalar diquark [Q, q] and 
axial vector diquark {Q, q} ground states are given in Table [0 

At the final step, we calculate the masses of heavy tetraquarks considered as the bound 
states of a heavy-light diquark and antidiquark. In this picture of heavy tetraquarks both 
scalar S (asymmetric in flavour [Qq] s=0 = [Qq]) and axial vector A (symmetric in flavour 
[Qq] s=1 = {Qq}) diquarks are considered. Therefore we get the following structure of the 
[Q9][0(f] ground (15) states (C is defined only for q = q')\ 



Two states with J PC = +A 



X(0++) = [Qq] s =o[Qq'}s=o 
X(0 ++ ') = [Qq] s=1 [Qq'] s=l 



Three states with J = 1: 



A^(l ++ ) = -^([Qq]s=i[Qq']s=o + [Qq}s= [Qq]s=i) 

X(l + -) = -^([Qq]s=o[Qq'}s=i ~ [Qq} s=1 [Qq'} s=0 ) 
X(l+-') = [Qq} s=1 [Qq'} s=1 
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TABLE II: Masses of charm diquark-antidiquark ground (15) states (in MeV). S and A denote 
scalar and axial vector diquarks. 



State 


Diquark 






Mass 




J PC 


content 




cqcq 


cscs 


cscq/cqcs 


0++ 


SS 




3812 


4051 


3922 


l+± 


{SA±SA)/V2 


3871 


4113 


3982 


0++ 


AA 




3852 


4110 


3967 


1+- 


AA 




3890 


4143 


4004 


2++ 


AA 




3968 


4209 


4080 


TABLE III: Thresholds for open charm decays and nearby hidden-charm thresholds. 


Channel 


Threshold (MeV) 


Channel 


Threshold (MeV) 


Channel 


Threshold (MeV) 


D°D° 


3729.4 


DfD- 


3936.2 


D°Df 


3832.9 


D+D- 


3738.8 


rj J '/ip 


4054.7 


D±DJ 


3837.7 


D°D*° 


3871.3 




4080.0 


D *o D ± 


3975.0 


pJ/4> 


3872.7 




4116.4 


D°D* ± 


3976.7 




3879.5 


D*+D* s - 


4223.8 




3988.6 




3879.6 






K*°J/ip 


3993.0 




4013.6 






D *o D *± 


4118.8 



• One state with J PC = 2 ++ : 

X(2++) = [Qq] s=1 [Qq'} s=1 . 

The orbitally excited (IP, ID . . .) states are constructed analogously. As we see a very rich 
spectrum of tetraquarks emerges. However the number of states in the considered diquark- 
antidiquark picture is significantly less than in the genuine four-quark approach. 

The diquark-antidiquark model of heavy tetraquarks predicts [6fl the existence of the 
flavour SU (3) nonet of states with hidden charm or beauty (Q = c,b): four tetraquarks 
( [Qq] [Qo~] > q = u,d) with neither open or hidden strangeness, which have electric charges 
or ±1 and isospin or 1; four tetraquarks ([Qs][Qg] and [Qg][Qs], q = u,d) with open 
strangeness (S = ±1), which have electric charges or ±1 and isospin |; one tetraquark 
([Qs][Qs]) with hidden strangeness and zero electric charge. Since in our model we ne- 
glect the mass difference of u and d quarks and electromagnetic interactions, corresponding 
tetraquarks will be degenerate in mass. A more detailed analysis predicts that such mass 
differences can be of a few MeV so that the isospin invariance is broken for the [Qg] [Qg] 
mass eigenstates and thus in their strong decays. The (non)observation of such states will 
be a crucial test of the tetraquark model. 

The calculated masses of the heavy tetraquark ground (IS) states and the corresponding 
open charm and bottom thresholds are given in Tables [ITW1 We find that all S"-wave 
tetraquarks with hidden bottom lie considerably below open bottom thresholds and thus they 
should be narrow states which can be observed experimentally. This prediction significantly 
differs from the molecular picture where bound B — B* states are expected to lie very close 
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TABLE IV: Masses of bottom diquark-antidiquark ground (IS) states (in MeV). S and A denote 
scalar and axial vector diquarks. 



State 
J PC 


Diquark 
content 




Mass 




bqbq 


bsbs 


bsbq/bqbs 


0++ 


SS 


10471 


10662 


10572 


l+± 


(SA ± SA)/y/2 


10492 


10682 


10593 


0++ 


AA 


10473 


10671 


10584 


1+- 


AA 


10494 


10686 


10599 


2++ 


AA 


10534 


10716 


10628 



TABLE V: Thresholds for open bottom decays. 



Channel 


Threshold (MeV) 


Channel 


Threshold (MeV) 


Channel 


Threshold (MeV) 


BB 


10558 


B+B- 


10739 


BB S 


10649 


BB* 


10604 


BfB** 


10786 


B*B S 


10695 


B*B* 


10650 


B*+B* s ~ 


10833 


B*B* 


10742 



(only few MeV below) to the corresponding thresholds. 

The situation in the hidden charm sector is considerably more complicated, since most of 
the tetraquark states are predicted to lie either above or only slightly below corresponding 
open charm thresholds. This difference is the consequence of the fact that the charm quark 
mass is substantially smaller than the bottom quark mass. As a result the binding energies 
in the charm sector are significantly smaller than those in the bottom sector. 

In Table fVTI we compare our results (EFG [4J) for the masses of the ground and excited 
charm diquark-antidiquark bound states with the predictions of Ref. p, Q, @, 0] and with 
the masses of the recently observed excited charmonium-like states [l|. We assume that the 
excitations occur only inside the diquark-antidiquark bound system. Possible excitations 
of diquarks are not considered. Our calculation of the heavy baryon masses supports such 
scheme In this table we give our predictions only for some of the masses of the or- 
bitally and radially excited states for which possible experimental candidates are available. 
The differences in some of the presented theoretical mass values can be attributed to the 
substantial distinctions in the used approaches. We describe the diquarks dynamically as 
quark-quark bound systems and calculate their masses and form factors, while in Ref. 0] 
they are treated only phenomenologically. Then we consider the tetraquark as purely the 
diquark-antidiquark bound system. In distinction Maini et al. consider a hyperfme inter- 
action between all quarks which, e.g., causes the splitting of 1 ++ and l + ~ states arising 
from the SA diquark-antidiquark compositions. From Table IVII we see that our dynamical 
calculation supports the assumption p that X(3872) can be the axial vector 1 ++ tetraquark 
state composed from the scalar and axial vector diquark and antidiquark in the relative IS" 
state. Recent Belle and BaBar results indicate the existence of a second X(3875) particle a 
few MeV above X(3872). This state could be naturally identified with the second neutral 
particle predicted by the tetraquark model [2j. On the other hand, in our model the lightest 
scalar ++ tetraquark is predicted to be above the open charm threshold DD and thus to 
be broad, while in the model [6| it lies few MeV below this threshold, and thus is predicted 
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TABLE VI: Comparison of theoretical predictions for the masses of the ground and excited charm 
diquark-antidiquark states cqcq (in MeV) and possible experimental candidates. 



orare 


Diquark 




Theory 


Experiment 


T PC 

j 


content 




Maiani et al. Maiani et al. (cscs) 


state 


mass 


IS 












n++ 

U 


qq 


^81 9 








1++ 


{SA + SA)/V2 


3871 


3872* 


f X(3872) f 
\ A"(3876) \ 


3871.4 ±0.6 
3875.4 ± 0.7+20 


i H — 
1 


{Dsi — D si) / y Z 


QS71 


O f 04 






n++ 

u 


A A 




oooz 






1 


A A 


oonn 

ooyU 


QQQO 

oooz 






2++ 


AA 


3968 


3952 


1 

y(3943) <^ 

1 


' 3943 ± 11 ± 13 
3914.3+H 

v. — o.o 














1— 


SS 


4244 


4330±70 


y(4260) 


f 4259 ± 8 + c 
1 4247 ± 12 + il 


1 

1 


( Q A Q A \ / /o 

[DJi — oA.) [y I 


A OQ A 

4ZS4 




y(4260) 


4283ljg±4 


1 

1 


A A 


/1 977 
4Z M 








1" 


AA 


4350 




y(4360) 


4361±9±9 


25 












l +± 


{SA±SA)/V2 


4431 




Z(4430) 


4433±4±1 


0++ 


AA 


4434 








1+- 


AA 


4461 


~ 4470 






2P 












1~ 


SS 


4666 




y(4660) 


4664±11±5 



' input 



to be narrow. Our 2 ++ state also lies higher than the one in Ref. [6|, thus making the inter- 
pretation of this state as y(3943) less probable especially if one averages the original Belle 
mass with the recent BaBar value wich is somewhat lower. 

The recent discovery of the y(4260), y(4360) and y(4660) indicates an excess of the 
expected charmonium 1 states [lj. The absence of open charm production is also in- 
consistent with a conventional cc explanation. Maini et al. [s| argue that y(4260) is the 
1 IP state of the charm-strange diquark-antidiquark tetraquark. We find that y(4260) 
cannot be interpreted in this way, since the mass of such ([cs]s=o[cs]s=o) tetraquark is found 
to be ~ 200 MeV higher. A more natural tetraquark interpretation could be the 1 IP 
state ([cq , ]s = o[cg]s=o) (SS) which mass is predicted in our model to be close to the mass of 
y(4260) (see Table IVT]) . Then the y(4260) would decay dominantly into DD pairs. The 
other possible interpretations of y(4260) are the 1 IP states of (SA — SA)/\/2 and AA 
tetraquarks which predicted masses have close values. These additional tetraquark states 
could be responsible for the mass difference of Y (4260) observed in different decay channels. 
As we see from Table I VI I the recently discovered resonances Y (4360) and Y (4660) in the 
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e + e~ — > 7t + 7t~iI)' cross section can be interpreted as the excited 1 IP (AA) and 2P (SS) 
tetraquark states, respectively. 

Very recently the Belle Collaboration reported observation of a relatively narrow en- 
hancement in the 7r + ?// invariant mass distribution in the B — > Kix + ip' decay [l|]. This new 
resonance, Z + (4430), is unique among other exotic meson candidates since it has a non-zero 
electric charge. Different theoretical interpretations were suggested [lj. Maiani et al. [q| 
give qualitative arguments that the Z + (4430) could be the first radial excitation (25) of a 
diquark-antidiquark X^(l + ~; IS) state (AA) with mass 3882 MeV. Our calculations indi- 
cate that the Z + (4430) can indeed be the l + ~ 2S [cu][cd] tetraquark state. It is the first 
radial excitation of the ground state (SA — SA)/y/2, which has the same mass as X(3872). 

In summary, we calculated the masses of heavy tetraquarks with hidden charm and bot- 
tom in the diquark-antidiquark picture. In contrast to previous phenomenological treatments 
we used the dynamical approach based on the relativistic quark model. Both diquark and 
tetraquark masses were obtained by numerical solution of the quasipotential equation with 
the corresponding relativistic potentials. The diquark structure was also taken into account 
with the help of the diquark-gluon form factor expressed in terms of diquark wave functions. 
It is important to emphasize that, in our analysis, we did not introduce any free adjustable 
parameters but used their fixed values from our previous considerations of heavy and light 
meson properties. It was found that the X(3872), F(4260), F(4360), Z(4433) and F(4660) 
exotic meson candidates can be tetraquark states with hidden charm. The ground states of 
bottom tetraquarks are predicted to have masses below the open bottom threshold and thus 
should be narrow. 
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